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ABSTRACT 

Flux-limited X-ray samples indicate that about half of rich galaxy clusters have cool cores. Why do 
only some clusters have cool cores while others do not? In this paper, cosmological N-body + Eulerian 
hydrodynamic simulations, including radiative cooling and heating, are used to address this question 
as we examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. These 
adaptive mesh refinement simulations produce both CC and NCC clusters in the same volume. They 
have a peak resolution of 15.6 /i -1 kpc within a (256 h^Mpc) 3 box. Our simulations suggest that 
there are important evolutionary differences between CC clusters and their NCC counterparts. Many 
of the numerical CC clusters accreted mass more slowly over time and grew enhanced cool cores 
via hierarchical mergers; when late major mergers occurred, the CC's survived the collisions. By 
contrast, NCC clusters experienced major mergers early in their evolution that destroyed embryonic 
cool cores and produced conditions that prevented CC re-formation. As a result, our simulations 
predict observationally testable distinctions in the properties of CC and NCC beyond the core regions 
in clusters. In particular, we find differences between CC versus NCC clusters in the shapes of X-ray 
surface brightness profiles, between the temperatures and hardness ratios beyond the cores, between 
the distribution of masses, and between their supercluster environs. It also appears that CC clusters 
are no closer to hydrostatic equilibrium than NCC clusters, an issue important for precision cosmology 
measurements . 

Subject headings: cool cores — galaxies: clusters: general — cosmology: theory — hydrodynamics — 
methods: numerical — intergalactic medium 



1. INTRODUCTION 

Clusters of galaxies with "cool cores" have peaked X- 
ray emission (i.e., excess above that extrapolated in- 
ward from a /3-modcl fit of the X-ray profile beyond the 
core) coincident with supergiant elliptical galaxies. They 
have central cooling times typically < 0.1i/ o _1 and cen- 
tral g as temperatures ~30-40% of the virial temperatures 
(e.g., Hkebe et alJll997tlLewis et al.ll2002UPeterson et ail 
2003). Cores of cool gas are found to be common in flux- 
limited samples, although a selection bias is likely present 
because of the strongly peaked X-ray emission in these 
clusters. Fro m a sample c o mpos ed of clusters detected 
with Einstein, IWhite et alj (|1997f ) found cool cores (CC) 
in «60% of their 207 cluster sample. Fro m a sample 
of 55 flux-limited RO SAT-observed clusters. iPeres et al.l 
((1991 found CCSs in over 70% of their galaxy clusters. 
More recently, IChen et al.l (|2007h identified 49% of their 
106 clusters as having cool cores in a flux-limited sam- 
ple, HIFLUGCS, based upon both ROSAT and ASCA 
observations. 

Why do some, but not all, galaxy clusters contain cool 
cores? To answer this question, we must explore the ori- 
gin and evolution of cool cores galaxy clusters. The earli- 
est and simplest model assumed clusters to be spherical, 
isolated systems where "cooling flows" formed; as radi- 
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ating gas loses pressure support, cooling gas flows in- 
wards to higher densit y values whic h further accelerates 
the cooling rate (e.g., lFabianl[2002f K But, the predicted 
end-products of this mass infall (e.g., star formation, HI, 
CO) have not been observed and the central tempera- 
tures indicate that the gas at the cores has only m oder- 
ately cooled (see review bv lDonahue fc Voitll2004l ). The 
current paradigm calls for heating to offset cooling, pos- 
sibly by AGNs via stro ng shocks (e.g.. iHeinz et al.lll998| ) 
or weak shocks (e.g., iFabian et al.l I2003T). or bv AG Ns 
+ thermal conduction (e.g.. iRuszkowski et al.ll2004f) . o r 
by AGNs + preheating (e.g.. iMcCarthv et alj|2007al ). 
or by Alfven waves induced by AGNs in t he inner core 
and c luster mergers in the outer cores (e.g., iFujita et al.l 
l200l . 

The simple cooling flow model did not incorporate the 
important effects of mergers and on-going mass accre- 
tion from the super cluster environment in which thes e 
clusters reside (e.g ., iMotl et al.ir2004t lPoole et al.l l2006l). 
iBurns et all (|1997h and iGomez et al.1 (|2002l) first exam- 
ined the result of cluster collisions in 2-D numerical 
simulations that involved two idealized spherical clus- 
ters with /3-model density profiles and central cool cores 
that collided together head-on. They found that the 
ram pressure from major mergers (i.e., subcluster to 
cluster mass ratios of »1 5% to 100%) tends to di srupt 
the cool cores. Similarly , iRitchie fc Thomasl (|2002l ) and 
iRicker &: Sarazinl (|2001h found disruptions of cool cores 
by major mergers between spherical clusters using 3-D 
simulations. These numerical models may suggest that 
the numbers of cool cores diminish as clusters grow via 
mergers (i.e., fewer cool cores in richer clusters at smaller 
z). 

More recently, we performed numerical simulations 
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of the formation and evolution of clusters in a cosmo- 
logical context using the adaptive mesh refinement N- 
body /hydr o code Enzo, aimed at further understanding 
cool cores (|Burns et al.ll2004 iMotl et al.ll2004f ). The gas 
in these clusters was evolved with radiative cooling but 
no heating. We found that cooling modifies not only 
the cores but also significantly alters t he cluster appear- 
ance out to the virial radius (see also lAkahori fc Masail 
2006). As new subcluster halos fall into a cluster, they 
gradually donate cool gas so that the cool cores grow 
over time. Most mergers are oblique with halos spiraling 
into the cluster centers and gently bequeathing cool gas 
to enhance the cores. Thus, in this model, cool cores 
themselves grow hierarchically via the merger /accretion 
process. This model predicts that even cool core clusters 
should possess a variety of substructures such as bul- 
let subclusters and cold fr onts, similar to those observed 
in Abell clusters (e.g. , lHallman fc Markevitchl 120041 : 
iMarkevitch et all 120021 : IMarkevitch fc V ikhlinin 120071 ). 
It also suggests cool cores may grow stronger (i.e., cooler 
and denser) as rich clusters increase in mass at recent 
epochs. 

The iMotl et all l|2004 ) simulations were limited by 
the baryonic physics that included only radiative cool- 
ing. This model suffers from the wel l-known cooling 
catastrophe (e.g., IWhite fc Reesl 119781 ) that results in 
an overproduction of cool cores and an inc rease in the 
baryon fraction (e.g., iKravtsov et al.|[2005l ). Nearly ev- 
ery dark matter potential well in this simulation was 
occupied by a halo of gas that had cooled signifi- 
cantly. Furthermore, these cool cores are "hard", gen- 
erally denser, colder, and with more distinct boundaries 
than are observed. The steep density contrast shel- 
ters the cores from ram pressure stripping, thus allow- 
ing them to survive and grow robustly during mergers 
with other halos. Clearly, a more realistic model of 
cool cores must involve added physical processes that 
"soften" the cores, thus making some susceptible to dis- 
ruptio n during mergers . Heating by star formation 
(e.g., lValdarnimll2006bD or by AGNs would potentially 
soften the cores. Additional sof tening effects may in- 
clude thermal conduction (e.g.. IZakamska fc Naravar] 
l200l iRuszkowski fc Beeelmanl [2002) along with mix- 
ing/heating from radio jet/lobe entrainment, and weak 
shocks and turbulent heating arising from halo mergers 
(e.g.. lBurnslll998t iBriiggen fc Kaiserll2002t iFuiita et al.l 
12004 IVoit"fc Donahudl2005l : iMathews et al 112006^ 

In this paper, we present Enzo cosmology simulations 
that include radiative cooling, star format ion (i.e., a 
mass sink for cold gas), an d heating (see also Mot l et all 
2005; H allman et al] |2006). Unlike previous simulations, 
our cooling + heating prescription has succeeded in pro- 
ducing both cool core and non-cool core clusters within 
the same computational volume. In addition to the some- 
what more realistic baryonic physics, these simulations 
have the advantage of bigger volumes and larger sam- 
ples o f clusters than in previous computatio nal simula- 
tions (|Motl et al J 12004 IKravtsov et al.ll2005l ). Thus, we 
have the dataset to examine evolutionary effects in these 
numerical clusters and can address the question in the 
title of this paper with good statistics. In Section 2, we 
describe the new simulations and the analysis of the nu- 
merical clusters. In Section 3, we compare the statistical 
properties of our numerical clusters with recent observed 



samples and show that the agreement is good. In Sec- 
tion 4, we describe new insights into the formation of 
cool core (CC) and non-cool core (NCC) clusters from 
our simulations. In Section 5, we describe the observa- 
tional consequences of evolutionary differences in CC and 
NCC clusters. Conclusions and a summary are presented 
in Section 6. 

2. NUMERICAL SIMULATIONS 

The simulations descri bed in this paper w ere performed 
with the Enzo 5 code (jO'Shea et al.ll2004h that couples 
an N-body algorithm for evolving the collisionless dark 
matter particles with an Eulerian hydrodynamics scheme 
(PPM) that utilizes adaptive mesh refinement. We adopt 
a ACDM cosmology with Sl b = 0.026, Q m = 0.3, fl A 
= 0.7, h = 0.7, n s =\, and erg = 0.9. The simulation 
was i nitialized at z=30 using the CDM transfer function 
from rEiscnstci n~fc Hul (|1999f l. A low-resolution simula- 
tion was first used to identify clusters in a volume of 256 
ft, -1 Mpc on a side using 128 3 dark matter particles and 
grid zones. High-resolution simulations were then per- 
formed that evolved the entire volume but adaptively re- 
fined 50 smaller regions (as separate simulations) around 
the largest clusters identified on the low resolution grid. 
Each of these 50 regions first is refined with two levels 
of static nested grids, each having a cell size half that of 
its parent grid (thus each region has spatial resolution 
4 times better than the parent grid). Within the nested 
grids, the dark matter particles have a mass resolution of 
9 x 10 9 h~ 1 M@. Then, within the nested static grids, we 
evolve the simulation with 5 additional levels of adaptive 
refinement, again with a factor of two increase in spatial 
resolution at each level. Cells are flagged for refinement 
based on the local baryonic and dark matter overden- 
sities, refining on thresholds of 8.0 times the minimum 
value at that level. The spatial resolution on the finest 
grid is 15.6ft. -1 kpc, adequate to resolve the cool core 
(rcc ~ lOO/i -1 kp c), but not to pr obe the details of its 
structure (see also lMotl et~aT1l2004) . 

Radiative cooling is calculated from a tabulated cool- 
ing curve der ived from a Raymond-S mith plasma emis- 
sion model (|Brickhouse et alj [[995 • assuming a con- 
stant metallicity of 0.3 relative to solar. The cooling 
curve is truncated below a temperature of 10 4 K. Every 
timestep, we calculate the energy radiated from each cell 
and remove that amount of energy from the cluster gas 
ijMotl et al.ll200l . 

As mentioned above, star formation provides one 
mechanism to soften cool cores by both transforming 
rapidly cooling gas into star particles (and, therefore, 
removing the cold gas) and by heating the surrounding 
gas with energy injected from supernovae. The star for- 
mation and heat ing that we used follow s the prescrip- 
ti on outlined by |Cen fc Ostrikerl (|1992l ) and described 
m iBurns et all (j2004f ). In brief, the code examines all 
grid cells at the finest refinement level above a spec- 
ified overdensity. The gas is converted to collisionless 
"star" particles when it is undergoing compression, rapid 
cooling, and the mass in the cell exceeds the JeanSs 
mass. The star formation rate is coupled to the local 
dynamical time and to a user-specified star formation 
efficiency. Once formed, the new star particle deposits 

5 http:/ /lea. ucsd.edu/portal/software/enzo 
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energy in the gas to simulate the instantaneous feedback 
from Type II supernovae. The strength of the supernova 
feedback is controlled by another efficiency parameter, 
e, which gives the thermal energy injected in proportion 
to the estimated rate of star formation for that particle 
(M s tar = M s tarltdyn, and e = eM star c 2 ). The most im- 
portant parameter in the star formation recipe was found 
to be the strength of thermal feedback from prompt su- 
pernovae. Through trial and error, we found that the 
value of e = 4.11 x 10~ 6 yields a reasonable fraction of 
bar yons locked in star particles at the current epoch (see 
also lBurns et al.ir2004h . This value for the feedback pa- 
rameter corresponds to (for a star formation rate of one 
solar mass per year) a supernovae rate of one per century 
with an average energy generation of 7 x 10 50 ergs per 
supernova. The chosen star formation parameters also 
produce both CC and NCC clusters in the same volume. 

In Figure [IJ we show a representative example of a 
cool core cluster at z = including images of the 
bolometric X-ray surface brightness, emission-weighted 
temperature, and th e distribution of star particles (see 
lHallman et al.l (|2006f ) for details on construction of syn- 
thetic X-ray and temperature images). These images il- 
lustrate the dynamic range in X-ray structures and tem- 
peratures typical in simulations of CC clusters, including 
the off-center infall of lower mass cool halos with leading 
bow shocks. The star particle image shows the distri- 
bution of sinks of cold gas and extended heating as new 
halos are accreted. 

The average total energy injection rate for the 10 most 
massive CC and 10 NCC clusters with a comparable mass 
distribution (at z = 0) within our computational vol- 
ume is w 5 x 10 43 ergs/sec. This is comparable to the 
X-ray luminosity for these clusters, thus our prescrip- 
tion produces an approximate balance between heating 
and cooling. This energy injection rate is also similar to 
the typical kinetic luminosities thought to powe r radio 
jets/l obes in central cluster radio sour ces {e.g.. iBurns 
fl99fllEilek fc Owenll2006l : lGentile et al.ll2007tlWlseet~aT 
12007ft . Thus, we view this energy injection scheme as a 
generalization of pre-heating of the core from a variety 
of sources, including AGN. The average energy injection 
rate is approximately the same over all epochs between 
z w 1 and z = 0. This feedback is also comparable for 
CC and NCC clusters of the same mass suggesting, as we 
show in Section 4, that something other than feedback, 
namely mergers, drive the evolution of these two cluster 
types. 

We constructed a catalog of all numerical clusters with 
-M200 > 10 14 M© from z=0 to z=2. (M200 is measured 
out to the radius, r2oo, where the density is 200 times 
the critical density and is ~ M V i r i a i which we will use 
interchangeably throughout this paper.) At z = 0, we 
have 94 clusters in the sample, but the entire catalog out 
to z — 2 contains 1522 clusters (many are the same clus- 
ter but at different epochs) giving us one of the largest 
samples of numerical rich clusters to date from a cosmo- 
logical simulation. These cluster simulations are publicly 
archived 6 . We have constructed a master table of the 
basic properties of these clusters including the average 
emission-weighted temperatures, virial and gas masses, 
f20o (~ the virial radius), baryon fractions, (5- model fit 

6 http://lca.ucsd.edu/data/sca 



parameters, CC or NCC designation, and other prop- 
erties. The ar c hive a nd this table will be presented in 
Hall man et al.l (|2007f ). These archived clusters form the 
basis of the analysis of numerical clusters presented in 
this paper. 

After visually inspecting all the temperature profiles 
for the z = clusters, we defined a cool core cluster to 
be one that has a >20% reduction in the central tem- 
perature compared to the surrounding region (where the 
slope of the temperature profile becomes negative; see 
Figure [7]) and this candidate baryonic cool core is within 
one zone («16 hr 1 kpc) of the dark matter density peak. 
This is a conservative definition that will yield the small- 
est number of cool cores, but we estimate that more 
liberal definitions will not increase the number by more 
than Ki 10%. With this strict classification, we find that 
16% of all the numerical clusters with M200 > 10 14 Af© 
have cool cores. This is lo w relative to the m ost recently 
observed fraction of 49% (|Chen et al.l 120071 ) that comes 
from a flux-limited sample (which may be biased some- 
what high by flux boosting from the cool cores). Several 
possible effects may be operating to reduce the fraction of 
numerical CC clusters. First, our baryon fraction (2.6%) 
for this simulation is now recognized as l ow relative to 
the re cent value from WMAP III (4.2%, ISpergel et all 
(2007)). A higher gas fraction could result in more ro- 
bust cool cores. Second, the power spectrum normal- 
ization (0%) may play a role in determining the number 
of CC clusters (our current value of cr 8 is larger than 
that inferred from WMAP III). Third, the numbers of 
cool cores and their survival during mergers appear to 
be a sensitive function of the heating/cooling prescrip- 
tion. Fourth, numerical resolution is likely a factor in 
the production of cool cores. 

3. STATISTICS OF THE X-RAY PROPERTIES OF 
CC AND NCC NUMERICAL CLUSTERS 

How well do the general characteristics of our numer- 
ical clusters match real galaxy clusters? This is an im- 
portant question to address before we propose a new for- 
mation scenario for cool core and non-cool core clusters 
based upon our numerical simulations. 

In the analysis that follows, we calculated the pro- 
jected average temperatures for our simula ted clusters 
as "sp ectroscopic-like" temperatures as in iRasia et al.l 
(2005). The weighting of temperature is different from 
the standard emission-weighted temperature, and has 
been shown to be more consistent with the value of 
the temperature which would be deduced from an X-ray 
spectral fit. The calculation performs 

500SL Jn 2 T a /r 3 / 2 dU' W 

where a=0.75 is the value fitted from FMazzott a et al.l 
(2004) which best approximates the value of the spec- 
troscopic temperature from X-ray fitting. In our case we 
have integrated this weighting in a cylinder with a radius 
of rsoo around the cluster center. 

In Figure [H the distributions of A/200 and T^qqsl 
(spectroscopic-like temperature inside r 500 ) for the over- 
all catalog of clusters, as well as for CC and NCC clusters 
within the catalog, are shown. We eliminated the cool 
core regions in calculating T500SL so as not to bias these 
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Fig. 1. — Representative example of a cool core cluster with 
M 2 oo = 5 X 10 14 Af Q and T virial = 4.2 keV at z = from the 
AMR simulation volume. The top image is projected synthetic 
X-ray surface brightness. The middle image is emission-weighted 
temperature (blue is T < 4 keV and yellow is T > 5 keV). The 
bottom image is a projection of the star particle density. The field 

eempM&^meMRfrements and to use a technique simi- 
lar to that applied to observations. As will be discussed 
further in Section 5, there are fewer high mass, high tem- 
perature CC clusters in comparison to NCC clusters in 
the catalog. At z = 0, the average mass of the CC clus- 
ters is 2.4 ± 1.4 x 10 14 M Q and that of the NCC clusters 
is 4.7 ±3.4 x 10 14 Af Q . 

In Figure we show examples of the most reliable 
statistics and basic relationships between variables that 
are typically calculated from X-ray observations. We 
compare these numerical data for z < 1 clusters (to ap- 
proximately match the range of redshifts for current ob- 
servations) in our catalog with those obtained from the 
recent statistically complete sample of clust ers observed 
with RO SAT and A SCA as reported by IChen et al.1 
(|2007f ). The numerical clusters were separated into CC 



and NCC using the criteria noted in Section 2. In the 
top panel of Figure 02 we plot the core radius versus the 
slope (/?) for a /3-model fit to the synthetic X-ray surface 
brightness profile (Sx oc [1 + (r/r c ) 2 ] 1 / 2-3/3 ). For the 
CC clusters, the cool cores were excluded from the fit 
(see Section 5.2 for details). This plot shows a separa- 
tion between CC and NCC clusters such that cool core 
clusters have smaller cores, r c , for a given 0. This sep- 
aration and the overall results fro m these /3-model fi ts 
match up very well with Figure 3 in IChen et al.l (|200l . 

The second panel shows the mean gas fraction mea- 
sured OUt to r 50 o (/ S as( r 50u) = (Pgas/ Ptotal) 5Qo) &S & 

function of T500SL ■ We attempted to correct our gas frac- 
tions by multiplying f gas by the ratio of f2f, measured by 
WMAP III to that which we used in these simulations 
(Section 2). This brings our gas fractions into somewhat 
better agreement but we emphasize th at they are still 
low relative to recent observat i ons (e.g., [M cCarth y et al.l 
I20n7btlvikhlinin et al.ll2nOfil3 : iSadat et al.ll2005r ). We do 
find good qualitative agreement in the shape and distri- 
bution of points in this figure relative to Figure 13 in 
IChen et al.1 (2007). There may be a slight tendency for 
reduced gas fr actions at lower t emperatures, as in obser - 
vations {e.g., iLin et all [20031 ; iMcCarthv et al.ll2007bD . 
but fg as is otherwise constant for T > 3 keV. There is a 
hint of a weak separation between CC and NCC clusters 
with cool cores having somewhat higher gas fractions for 
a given temperature as we will discuss in Section 4. 

The third panel presents a plot of gas mass out to 
?"500 against T^qqsl- There is a strong scaling relation 
with comparable power-law slopes for each type of clus- 
ter (measured slope index of 1.61 ± 0.04 for CC clusters 
and 1.69 ± 0.01 for NCC clusters). This scaling rela- 
tion is qua li tative ly similar to that in Figure 11 from 
IChen et al.1 (|2007l ). although the slope is a bit steeper 
than that obs erved and expected for self-similar behavior 
(|Kaiserlll986D {i.e., M oc T 15 ). iVikhlinin et all (|2006bf ) 
find a flatter slope for the M-T relation for their sample 
of 13 cool core clusters in comparison to other authors 
who analyzed mixed samples with CC and NCC clusters 
(and different techniques for measuring mass and tem- 
perature). 

In a separate paper (jJeltema et al.ll2007h . we also show 
that the distribution of X-ray substructure within these 
clusters, as measured using power ratios, agrees with that 
observed from X-ray observations of nearby rich galaxy 
clusters. 

Overall, within the noted limitations of these simula- 
tions, the average properties and the relationships be- 
tween basic variables for the numerical clusters agree 
fairly well with X-ray observations. 

4. THE FORMATION OF CC AND NCC CLUSTERS 

With a relatively large sample of numerical clusters, we 
are able to explore the question of why some clusters have 
cool cores but others do not. These simulations indicate 
that the evolution of CC clusters has followed a different 
history in terms of accretion of mass from the cosmic web 
in comparison to NCC clusters. Although each cluster 
has its own unique rich and complex evolutionary path 
that depends upon its initial mass and the density of 
surrounding halos, there are some clear general trends 
within this larger dispersion that we see which appear to 
separate CC from NCC clusters. 
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Fig. 2. — Distribution of the cumulative fractions of total cluster 
masses and T500SL (spectroscopic-like temperatures within rsoo) 
for all clusters (solid) and for each cluster type in the catalog with 
M20Q > 10 14 Mq and z < 1 (to approximately match present range 
of X-ray cluster observations). 

In Figure [?J we track the evolution of the median 
changes in cluster mass (M200) and the median central 
temperatures for the 10 most massive CC and NCC clus- 
ters from z = back to between z = 1 and z — 2 (each 
cluster is tracked back only to a time determined by our 
mass cutoff of A/ 2 oo > 10 14 M Q ). Although we use the 
ten highest masses in each sample for better statistics, 
we note that the same qualitative trends as shown in 
Figure [H but with larger dispersion, are present for sam- 
ples of CC and NCC clusters selected to have comparable 
mass distributions. We illustrate this evolution back to 
just z w 1.5 because there are only a few CC clusters 
at z > 1.5 above our 10 14 M Q mass limit (more clusters 
grow above this mass over time). 

The two cluster types show different histories in their 
median mass accretion rates at early times. At z = 1.5, 
NCC clusters experience a median «75% change in mass 
per Gyr, albeit with a not unexpected large dispersion 
due to the wide range of merger states. The CC clusters 
have a median change of ~30% in mass per Gyr with a 
smaller dispersion but also for fewer clusters with M200 > 
10 14 M Q . Using a K-S test, we find that for 1 < z < 1.5, 
the distributions of mass change for NCC and CC clusters 
differ at the 95% level. 

The NCC clusters demonstrate a trend of experienc- 
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Fig. 3. — Statistical properties of numerical galaxy clusters with 
M200 > 1O 14 M and z < 1. Blue are CC and red are NCC clus- 
ters. Top: /3-model fits to Sx profiles of individual clusters pro- 
vide a measure of the core radius (r core ) and slope (/3). Middle: 
Gas fraction versus T500SL (spectroscopic-like temperatures within 
rsoo)- Bottom: Gas mass out to rsoo is plotted against TsooS-L- 
Power-law fits, performed separately for CC and NCC clusters, are 
shown. 

ing major mergers early in their histories up to z w 0.5, 
which destroy any initial cool cores, then they settle down 
to a more quiescent state thereafter. Here we define a 
"major" merger as one that has the potential for disrupt- 
ing a nascent cool core, usually accreting >50% of the 
cluster's previous mass over a timescale of «1 Gyr. CC 
clusters, on the other hand, avoid mergers with high frac- 



6 



tional mass changes early in their histories and instead 
grow slowly such that the cool cores increase in mass and 
stability. As shown in Figure^] CC clusters after z w 0.5 
have a relatively constant rate of accretion continuing to 
the present, similar to NCC clusters. 

The central temperature plot in Figure [4] demonstrates 
that similar starting conditions can result in either CC or 
NCC clusters. At early epochs, the dispersion in central 
temperatures is large and the distributions are statisti- 
cally indistinguishable between what will become CC and 
NCC clusters at z — 0. This contrasts to the significant 
difference in central temperatures between CC and NCC 
clusters for z < 0.5. Thus, the early merger history pri- 
marily determines the eventual cluster configuration at 
the present epoch. 

In Figures O and [51 we show examples of the evolution 
of NCC and CC clusters, respectively, which well rep- 
resent the general scenarios for how these clusters form. 
Our simulations indicate that lower mass clusters with 
T < 2 keV form cool cores early in their history when 
initial conditions produce central densities and tempera- 
tures that allow the gas to radiatively cool. This suggests 
that many (most) lower mass clusters should have cool 
cores which is consistent with the data in Figure [5] in the 
next section and the observations of poor clu sters com- 
posed of early-type galaxies (see e.g . , review bylMulchaeyl 
' 20041: IZabludoff kz Mulchaevlll99St iPonman et alJl2003t 
Chen et all 120071) . According to our simulations, early 
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mergers cause the fates of NCC and CC clusters to di- 
verge. 

As shown by the representative example in Figure 
NCC clusters often undergo major mergers early in their 
history. This cluster had begun to develop a cool core 
at z = 1.25 (see panel A). However, the cluster expe- 
rienced a major merger (mass increased by ~100%) at 
2 w 1 and the cool core was greatly diminished (panel 
B). Smaller mass halos with cool cores continue to be 
accreted by this cluster but these CC's are ram pressure 
stripped/disrupted usually within a single core passage 
(see panels C and D). By z ~ 0.65, there is no evidence 
of a CC associated with the cluster dark matter density 
peak at this or later times. Early mergers destroy the 
cool cores in NCC clusters, leaving behind hotter, ther- 
malized, moderately dense cores where the cooling time is 
everywhere above the Hubble time. As shown in the next 
Section and in the Appendix, the NCC cluster gas has 
become mostly relaxed within the gravitational potential 
well (with minor perturbations from small infalling halos) 
with a surface brightness profile well represented by a (3- 
model. Subsequently, cool halos infalling into these NCC 
clusters do not survive passage through the central parts 
of the clusters nor do the central conditions allow cool 
cores to re-establish. NCC clusters continue to experi- 
ence minor mergers as they now slowly evolve (typically, 
mass increases only as 10% over Gyr time frames after 
z > 0.5 from multiple mergers for the NCC as shown in 
Figured]). We suggest that such an early major merger 
produced the characteristics observed today for the NCC 
Coma cluster whose com plex properties m ay be the re- 
sult of previous mergers (|Burns et al.lll994 h 

On the other hand, Figure [6] suggests that CC clus- 
ters evolve differently This CC cluster had no signifi- 
cant change in mass until z = 0.75 and its only major 
merger did not occur until z — 0.3. Figure [6] shows the 
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Fig. 4. — The median time evolution of the 10 most massive 
NCC clusters and the 10 most massive CC clusters. Dispersions 
in the distributions are also shown. Top: The median fractional 
mass change per Gyr as a function of lookback time and redshift. 
Bottom: Evolution of the median core temperatures (normalized 
by the virial temperature). The NCC clusters accrete significantly 
more mass than the CC clusters until z 0.5, signifying more early 
major mergers than for the CC clusters. By z = 0.25, the core 
temperatures for the NCC clusters are about 3 times hotter than 
the CC clusters; the cool cores are well established and becoming 
more robust (slightly cooler and denser) throughout subsequent 
minor mergers. j i i . i 

temperature and central cool core as tne merger is pro- 
gressing (at z = 0.3) (panel A). The next snapshot at 
z = 0.2 shows the CC somewhat diminished but still 
present. The central temperature moved slightly up- 
ward but quickly readjusted downward as the CC eas- 
ily survives the shock heating and ram pressure from the 
merger. In contrast to the simple cooling flow model, CC 
clusters may be no closer to hydrostatic equilibrium than 
NCC clusters with the equivalent mass (see also Section 
5.5). This hierarchical formation model for CC clusters 
makes clear predictions of substructure and average clus- 
ter characteristics beyond the core that are testable with 
X-ray data. 

Figure [7] shows the radial profiles of the baryon frac- 
tion and the temperature for representative examples of 
numerical CC and NCC clusters in our sample. Even 
outside of the cool core (ss 0.05r2oo ~ lOO/i -1 kpc, the 
first vertical dotted line), there is an excess of baryons 
relative to NCC clusters out to w 0.3r2oo (second vertical 
dotted line). (We note that the dark matter density pro- 
files are comparable for the CC and NCC clusters.) Such 
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Fig. 5. — Four snapshots of the history of an NCC cluster with final (. 
vertical lines in the bottom panel correspond to different epochs of the 

an extended "transition region" could be created, in part, 
by gas "sloshing" in the cluster p otential well following re- 
peated mergers as proposed by M arkevitch fc Vik hlinin 
(T2007h . 

Figure [7] also shows the temperature differences be- 
tween each cluster type. The NCC cluster demonstrates 
the universa l temp erature profile that we described in 
lLoken et al.l ((2002). The temperature profile of the CC 
cluster rises steeply to « 0.05r 2 oo and then it has a 



! = 0) values of M 2 oo = 8 X 10 14 Af Q and T virial = 5.5 keV. Dashed 
temperature images. 

prolonged stretch of near-constant temperature, again 
within the region « 0.05r 2 oo to « 0.3r2oo- This CC clus- 
ter profile is less compact than that found for recent SPH 
simulations (Valdarnini 20 06bf) and agrees well with ob- 
serva tions (see e.g.. I Vikhlinin et all 2006b: Baldi et all 
12001 . 

Using these two trends, we define three components 
to a cool core cluster: cool core, "transition region", 
and outer region. The transition region is differenti- 
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Fig. 6. — Two snapshots of the history of a CC cluster with M200 

ated by the excess of baryons outside the core and rel- 
atively flat temperature profile (and low entropy) com- 
pared to NCC clusters. We have (subjectively) chosen 
the limits ~ 0.05r2oo to « 0.3r2oo for this transition 
region; w 0.05r2oo is the traditional edge of the cool 
core where the slope of temperature changes dramati- 
cally, and ~ 0.3r2oo is an average location where the 
baryon fraction of each type of cluster converges and the 
temperature begins to decrease. Most current cluster X- 
ray observations also measure Sx accurately out to only 
~ 0.3r2oo, meaning that most observations measure pre- 
dominantly the transit ion region in CC clus ters (as will 
be discussed further in lGantner et al.l (|2007f )). 

These simulations predict a very different set of cluster 
characteristics from those expected in the simple, non- 
evolving cooling flow model or from cooling-only simula- 
tions. Since CC and NCC clusters have experienced dif- 
ferent magnitudes and epochs of mergers, there should 
be observational signatures remaining from the merg- 
ers. In particular, we will show in the next section that 
the fraction of clusters with cool cores is expected to be 
strong function of the virial mass, in good agreement 
with recent observations. We then will show that sin- 
gle /3-models systematically overestimate (or bias) the 
densities and masses beyond the cores in CC clusters. 
Furthermore, our simulations predict that more cool gas 
should be found beyond the cores in CC clusters in com- 



: 4.4 X 1O 14 M and T virial = 3.7 keV at 2 = 0. 

parison to NCC clusters. Finally, this scenario forecasts 
that CC rich clusters should be found in denser super- 
cluster environments at the present epoch. 

5. CONSEQUENCES OF EVOLUTIONARY 
DIFFERENCES IN CC AND NCC CLUSTERS 

In this Section, we explore the differences in the prop- 
erties of CC and NCC clusters based upon the results of 
our numerical simulations. The simulations predict sub- 
stantial differences in the characteristics of these clusters 
beyond the cores. These predictions can be tested with 
data from current and planned X-ray telescopes. 

5.1. Masses and Fractions of CC Clusters 

An intriguing new result is that the fraction of clus- 
ters with cool cores is a strong function of cluster gas 
mass as shown in Figure [5] We display gas masses here 
instead of total cluster masses to allow a direct com- 
parison with t he ob servat ion-derived data p resented by 
lO'Hara et ail (|2006f) and IChen et all (|2007l ). About a 
quarter of simulated clusters with M gas « 5 x 10 12 Af© 
have cool cores whereas no high mass numerical clusters 
{M gas > 4 x 10 13 M Q ) have cool cores. As a corollary to 
this result, we find that the mean total mass for the 10 
most massive CC clusters at z = is 2.4 ± 1.4 x 10 14 M Q 
whereas the mean total mass for the 10 most massive 
NCC clusters is 11.3 ±4.0 x 10 14 M o . These results sup- 
port the idea that cool cores are destroyed via multiple 
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Fig. 7. — Plots of the gas fraction and temperature profiles for 
representative CC and NCC numerical clusters. The vertical lines 
designate the location of the newly defined CC cluster transition 
region (0.05r2oo to 0.3r2oo)- This transition region contains an 
excess of baryons compared to a typical NCC cluster and a rela- 
tively fiat temperature profile, indicating different gas properties 
than either the core or the outer region. 

major mergers and the probability of cool core disruption 
increases as clusters grow to the size of the Coma cluster. 

A similar result can be seen from observa tional samples 
of clus ters usi ng ROSAT data co mpiled by lO'Hara et alj 
(|2006ft and bv lChen etai] (|2007t) . Both sam ples have a 
somew hat common ancestry from the work of lEdge et al.l 
(1990) with the samples consisting of nearby (0.01 < 
z < 0.1), moderate X-ray lumino sity clusters. I n Fig- 
ure p we have plotte d data from lO'Hara et alj (|2006f) 
and lChen et alJ (|2007l ) overlaid onto those for our sam- 
ple of numerical clusters. Although the absolute values 
of the fractions differ between the observed and numer- 
ical samples as discussed in Section 2, the general trend 
of decreasing fraction of CC clusters with mass is present 
for both observations and simulations. 

This finding is contrary to the expectations of the sim- 
ple non-evolving cooling flow model where the number of 
cool cores should increase with cluster mass (as the cen- 
tral gas density increases). Although obs ervational se- 
lection effe cts are a pos s ible co ncern in the lO'Hara et al.l 
(|2006T ) and I Chen et all (|2007l ) catalogs, we believe that 
this newly discovered trend contains important insights 
into the formation of CC versus NCC clusters. 



5.2. Evolution in the Fraction of Cool Cores? 

Recently, IVikhlinin et alj (|2006al ) reported that the 
number of observed CC clusters declines dramatically to 
15% with redshift beyond z sw 0.5 (versus 65% for their 
nearby cluster comparison sample). Because of the lim- 
ited spatial resolution with Chandra at these distances, 
they use the central slope or "cuspiness of the surface 
brightness" to distinguish between CC and NCC clusters. 

In Figure [HI we show the fraction within the co- moving 
volume of our numerical cool core clusters as a function 
of redshift for all clusters with M > 10 14 M Q out toz~ 1. 
The error bars in each bin reflect the y/~N uncertainties 
due to the number counts. Within these errors, the frac- 
tion of CC clusters is not a strong function of redshift 
out to z w 1 (15-20%). For z > 1, the fraction drops to 
1=3 10% but the dispersion is large because there are only a 
few CC clusters with M > 1O 14 M0 at these early epochs. 

The flat distribution of CC fraction within the range 
< z < 1 is not inconsistent with the evolutionary for- 
mation scenario described in Section 4 and shown in Fig- 
ure [J] for several reasons. First, Figure 2] reveals that the 
greatest disparity in mass change between CC and NCC 
clusters occurs for z > 0.75. That is, most of the growth 
in NCC clusters via mergers occurs at the expense of CC 
clusters at earlier epochs. For z < 0.75, CC and NCC 
clusters grow at comparable rates. Second, although CC 
clusters continue to be lost via mergers for z < 0.75, this 
is counterbalanced by the fact that the numbers of CC 
clusters above the mass cutoff of 1O 14 M0 continue to in- 
crease because CC clusters also grow via accretion. Thus, 
the rate at which CC clusters are destroyed is approxi- 
mately equal to the rate at which new clusters are added 
above our 1O 14 M mass limit. This produces the effect 
of no apparent evolution in the fraction of CC clusters 
for z < 1 in a mass-limite d sample. 

In order to explain the IVikhlinin et alj (|2006af ) result 
within the context of our simulations, their sample would 
have to have substantial selection effects possibly driven 
by the unique choice of cool cores based upon the slope 
of the X-ray surface brightness profile and resolution ef- 
fects. Alternatively, some form of time-dependent baryon 
physics (e.g., higher f eedback rates a t earlie r epochs as 
recently described by lEastman et alj (|2007f) for AGNs) 
not incorporated into the present simulations could po- 
tentially boost the fraction of CC clusters seen in the 
local Universe compared to earlier epochs. 

5.3. Surface Brightness Profiles 

Turning next to the large-scale X-ray surface bright- 
ness profiles (Sx), we fit /^-models in two different ways 
for two different subsamples of the numerical clusters to 
examine potential differences in the shapes and core radii 
of CC versus NCC clusters. 

We began by producing average Sx profiles from pro- 
jected X-ray images along a single axis for all the clus- 
ters at z — and with M > 1O 14 M0, separated into CC 
(10 clusters) and NCC (78 clusters) categories (several 
z = clusters were not used because of contamination 
by multiple clusters) . The flux for individual profiles was 
first normalized by h'hoo 3 ^ 2 (from the Mass- Temper ature 
scaling relationship as in iFino R'ucnov etHl (|200lh ) be- 
fore averaging them together. We then fit a /3-model to 
each of these average profiles as would be done for obser- 
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Fig. 8. — Fraction of all galaxy clusters with cool cores as a 
funct ion of mass. The s imulated data is from 1522 numerical clus- 
ters. [Chen et al.l ||2007|1 is bas ed upon their table of 106 observed 
clusters. IQ'HaraeTlrt] ||2006D is from their table of 45 observed 
clusters. Mgas^QQ was chosen as the measurement common across 
all three data sets. The differences in absolute fraction of CC clus- 
ters between the samples may have a variety of causes (e.g., lumi- 
nosity boosting by cool cores in flux-limited observational samples; 
numerical clusters may be affected by resolution effects, selection 
of cosmological parameters, or feedback prescription). 
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Fig. 9. — Fraction of cool co res as a function of red shift for 
numerical clusters. In contrast to Vikhlinin et al. (2006a), we find 
no change in the fraction of cool cores with redshift. 

vations. That is, we fit the profiles out to 0.3r2oo (roughly- 
corresponding to 0.5r 50 o ~ 0.5 Mpc) which is the typical 
limit to the observed surface brightness in most X-ray ex- 
p osures with cur r ent in struments (as we discuss further 
in lGantner et a l. (2007)). For the average CC profile, we 
excluded the cool core in making the /?-model fit. The 
result is shown in Figure [TOl 

This figure indicates that the profiles (beyond the cool 
core) for CC clusters are distinguished from NCC clus- 
ters in several important ways. First, the parameters 
for the (3- models are different. For the average CC pro- 
file, r c = 0.05 ± 0.09 r 200 and f3 = 0.66 ± 0.12, whereas 
for the NCC average S x profile, r c = 0.12 ± 0.02 r 200 
and j3 = 0.66 ± 0.07. That is, the cluster core radii are 
much smaller for CC clusters, as also shown in Figure [H 
There is also considerably more scatter in the fit for the 
CC average profile (consistent with more variation be- 



tween individual profiles) than for the NCC Sx profile. 
Second, the shape of the two Sx profiles are different 
within the transition region where the slope of the NCC 
is generally flatter than the CC cluster, as would be ex- 
pected from Figure [7] There are similar slope differences 
in the S x profiles between the NCC (Abell 401) and CC 
(Abell 85) clusters computed from deep Ch andra obser- 
vations reported bv lVikhlinin et aT] (|2006bf ). Third, the 
/3-model is a better fit to the average NCC cluster profile 
than to the CC profile. In particular, the /3-model fit to 
the regions that would be typically observed by current 
satellites (i.e., the transition region) for CC clusters sig- 
nificantly overshoots the actual flux in the outer parts of 
clusters. At r2oo, the /3-model overestimates the flux of 
the average CC profile by a factor of 3.8. This will result 
in a serious bias of cl uster gas masses as we discussed in 
iHallman et aT] (|2006l) . 

We also fit /3-models to individual profiles for all nu- 
merical CC and NCC clusters (from a single projection) 
in our master database with M > 5 x 1O 14 M0 and z < 2. 
In this case, we fit models out to r = r§oo . We did not use 
the inner portion of the profiles dominated by the cool 
core (determined by the point where the slope of the tem- 
perature profile becomes negative) in making fits to CC 
clusters. We then calculated the reduced x 2 goodness- 
of-fit values (including extrapolations of the fits out to 
^200) as compared with the numerical X-ray profiles. A 
histogram of those reduced \ 2 values for CC and NCC 
clusters is shown in Figure [11] As also indicated in Fig- 
ure [TOl the /3-models fit the NCC clusters much better 
than the CC clusters. About 88% of the NCC clusters 
have x 2 < 1, whereas about one-third of the CC clusters 
have x 2 > 1- Once again, this is caused by the slope 
changes from the transition region to the outer core in 
the CC clusters which is not fit well by a single /3-model. 

As we show in the Appendix, good /3-model fits to 
Sx suggest that a nonisothermal gas in the intraclus- 
ter medium (ICM) is in approximate (but not necessar- 
ily perfect, see Section 5.5) equilibrium with the gravita- 
tional potential well of the cluster (under the assumption 
that the ICM gas is polytropic). So, the above results 
suggest that the ICM in NCC clusters is approximated 
by a gas with a balance between heating and cooling that 
is quasi-relaxed in an NFW-like dark matter potential. 
This is consistent with Figure 2] which shows that NCC 
clusters underwent major mergers early in their history 
but have subsequently settled into a quasi-equilibrium 
state with only minor on-going mergers. This contrasts 
with CC clusters that have temperature and density pro- 
files inconsistent with a simple adiabatic gas well beyond 
t he cool core. 

Vikhlinin et al. (2006b) recently analyzed the gas and 
total mass profiles for 13 nearby, "relaxed" clusters with 
temperatures between 0.7 and 9 keV using data from 
Chandra. All of the clusters in this sample have cool 
cores. They examined the surface brightness profiles out 
to at least r5oo and concluded that the profiles are "not 
described well by a beta model". When single (3— models 
are fit to the inner portions of the clusters (but exclud- 
ing the cool cores) , they poorly extrapolate the gas den- 
sity and mass profiles in the outer parts of the CC clus- 
ters as we found for th e numerical clusters. In addition, 
Vikhl hiTn et alj (|2006b[ ) observed that the cooler regions 
in low temperature clusters are confined to a smaller frac- 
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Fig. 10. — Average synthetic X-ray surface brightness profiles 
for CC and NCC clusters with z = and M > 10 14 M Q (10 CC 
averaged together and 78 NCC averaged together). Error bars 
are errors on the mean determined from variations within each 
bin. Dashed curves are the best fit /3-models within the "transition 
region" (i.e., between the vertical lines from 0.05r200 to 0.3r20o) 
for the CC average cluster profile and the best fit including all 
the points out to the right-most vertical line for the NCC cluster 
profile. The NCC profile was arbitrarily shifted downward by a 
factor of 10 to better distinguish it from the CC profile. 




0.4 E 1 

2 4 6 8 10 

Reduced Chi Squared 



Fig. 11. — Histograms of the reduced \ 2 values for /3-model fits 
for CC (solid) and NCC (dashed) clusters. Each numerical cluster 
(M200 > 5 X 1O 14 M0, 2 < 2) is fit to rsoo and extrapolated out to 
7"200, and a x 2 goodness-of-fit is calculated for the entire profile. 

tion of the virial radius than in the hotter CC clusters. 
This is consistent with tr ends found in our nu merical 
simulations as described bv lHallman et all (|2006f ) and by 
lAkahori fc Masail (|2006f) . Finally, Vikhlinin et al. note 
that low temperature (mass) clusters (T < 2.5 keV) with 
cool cores have a bigger ratio of central to virial temper- 
ature than do the clusters with larger T V i r i a i. We plan 
to explore the origin of this effect with new higher reso- 
lution simulations with more sophisticated heating pre- 
scriptions. 

5.4. Temperatures and Hardness Ratios Beyond the 
Cluster Cores 

The temperature profiles for the two clusters shown in 
Figure [7] indicate that the temperature distributions for 



CC and NCC clusters are substantially different out to 
w 0.3r2oo- In particular, the broad "transition region" 
has both a cooler and flatter distribution of temperature 
outside the cool core for the CC cluster in comparison to 
the NCC cluster. How general is this result for the total 
sample of numerical clusters? To address this question, 
we made emission-weighted temperature images of the 
clusters in our numerical cluster catalog with M200 > 
5 x 10 14 M Q and < z < 0.5. From these images, we 
produced the histogram of temperatures (normalized by 
T v iriai), excluding cool cores, shown in Figure [T2l The 
distribution of temperatures beyond the cool cores for 
CC clusters is significantly different from NCC clusters 
with a broad tail toward lower temperatures. CC clusters 
have ^40% more gas with T ew < 0.3T„j r j o ( beyond the 
cores than NCC clusters. 

We predict that this signature will be apparent in hard- 
ness ratio maps that are commonly made from X-ray ob- 
servations. As shown in Figure [TBI for four typical cases 
drawn randomly from our simulations, the hard-to-soft 
band ratios (2-8 keV/0.5-2 keV) do a good job of illus- 
trating the abundance of cooler gas beyond the cores in 
CC clusters. Figure ITU shows the cumulative fraction of 
pixels below a given hardness ratio for all four clusters in 
Figure [T3l For the two CC clusters, we have excised the 
cool cores (< 0.05r2oo) so as not to bias the results with 
gas already known to be cooler than its NCC counter- 
part. As expected, the two CC clusters have a majority 
of pixels with values < 1 and therefore are cooler in the 
transition region than the NCC clusters. The NCC clus- 
ters are both centered approximately at hardness ratio 
ss 1, hence the gas in these clusters is roughly at the 
virial temperature o f the clusters. 

We shall show in iGantner et al.l (|2007t ) that there is 
very good agreement in the predicted hardness ratios 
from our simulations with X-ray observations of clusters 
from the Chandra and ROSAT archives. 

5.5. The Supercluster Environments of CC and NCC 

Clusters 

The mass and temperature evolution plots in Figure 
H] indicate that NCC clusters underwent major mergers 
early in their history in contrast to the milder accretion 
over time for CC clusters. This may also suggest that 
the larger scale environments in which these two types 
of clusters live are different since accretion of halos and 
diffuse material must come from the cosmic web. It is 
possible that NCC clusters began their lives in higher 
overdensity regions which then accele rated the growth of 
these clusters v ia mergers (see e.g., iMo fc Whitdll996l : 
iGao et al.|[2005] ). 

In an effort to explore the possible influence of the 
differences in the supercluster environments for CC and 
NCC clusters, we calculated the real space densities 
(which are expected to correlate with accretion rates) 
of all halos with M 20 o > 1O 13 M (« mass resolution of 
simulations) within a radius of 5r2oo of numerical rich 
clusters with virial masses 1 — 6 x 10 1 -M© (approximate 
mass range of the CC clusters as shown in Figure [5]). 
We calculated these densities for a series of redshifts be- 
tween and 1.5, and separated clusters by CC and NCC 
according to our definitions in Section 2. 

The ratio of supercluster densities for CC to NCC clus- 
ters as a function of redshift is shown in Figure [151 At 
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Fig. 12. — Histograms of the ratio of emission- weighted temper- 
atures (T ew ) from r = 180/i -1 kpc to r2oo divided by T viria i from 
temperature images for CC and NCC clusters. Clusters from the 
numerical catalog are included with M200 > 5 X lO 14 A/0. Red are 
NCC and blue are CC clusters. 

early epochs (z > 1), the average supercluster density 
is somewhat higher around NCC clusters than CC clus- 
ters. One might expect clusters that are experiencing 
major bouts of accretion of subclusters that result in 
the destruction of embryonic cool cores to be surrounded 
by a higher density of halos. At times corresponding 
to 0.7 < z < 1, there may a slight underdensity of halos 
around NCC clusters in comparison to CC clusters as one 
might expect if the NCC clusters suffered a large amount 
of mass accretion effectively "clearing out" its nearby 
neighborhood. Interestingly, at late times (z < 0.3), the 
density of halos around CC clusters is «30% greater than 
for NCC clusters. However, unlike earlier epochs where 
the mass ratio of the main cluster to the average neigh- 
boring cluster is often « a few, this average mass ratio of 
the rich cluster to the halos at z < 0.3 is much larger as 
the main cluster has grown considerably over the past 10 
Gyr. This means that there are many small subclusters 
falling into the CC clusters but their relative impact is 
small compared to that for NCC clusters at earlier epochs 
(see Figure [7]). However, the above trends are weak at 
best and there is a large dispersion in supercluster den- 
sities between individual clusters. 

It is interesting to note that this possible trend of over- 
abundance of halos around numerical CC cl usters at the 
presen t epoch is also found in Abell clusters. lLoken et alj 
(1999) constructed a volume-limited sample of z < 0.1 
Abell clusters that was estimated to be 98% complete in 
an effort to investigate their supercluster environs. They 
separated clusters into CC and NCC. They found that 
CC Abell clusters have twice the density of neighboring 
clusters as do NCC clusters out to a radius of 4 3 h -1 
Mpc. We attempted to mimic the lLoken et al.1 (|1999h 
analysis by recomputing halo densities out to radii of 43 
h^ 1 Mpc and included neighbor halos in the calculation 
only if they had masses > 10 14 M Q (i.e., Abell-like clus- 



ters). We found that the density of neighboring halos is 
«40% greater for CC than NCC clusters, slightly larger 
than in Figure [T51 but still considerably less than what 
lLoken et al.1 (|1999f ) propose for Abell clusters. 

5.6. Deviations from Hydrostatic Equilibrium 

Galaxy clusters are potentially powerful tools for 
precision cosmology. Accurate cluster mass determi- 
nations and gas fractions, along with cluster abun- 
dance counts, can provide key constraints on the dark 
energ y parameter, w, as wel l as fib, ^m, and as 
(e.g., IWang fe Steinhardtl[l99"8t lHaiman et al.il200lD . In 
recent efforts, samples of CC clusters are being used ex- 
clusively because the y are believed to b e dynamically re- 
laxed. For example. I Allen et al.l (|2007f ) selected a sam- 
ple of 42 hot, X-ray luminous clusters with 0.05 < z < 
1.1, all of which have short central cooling times (< 
a few 10 9 yrs), to constrain cosmol ogical parameters 
from f aas . Previo us simulations (e.g., iRasia et "al"1l2006l ; 
iNagai et al.ll2007l and references therein) have called into 
question hydrostatic equilibrium for clusters. But, are 
CC clusters really more dynamically relaxed than NCC 
clusters? Figures [4] and [15] seem to call this assumption 
into question. 

To explore this further, we calculated the devia- 
tions from hydrostatic equilibrium for all the clusters in 
our numerical arch ive with M 2 oo > 10 14 M© (see also 
Ueltema et al.1 [20071 ) . We did this by calculating the es- 
timated mass of clusters from the gradients in the tem- 
peratures and gas densities in the usual way assuming 
hydrostatic equilibrium. In this case, we have used the 
spherically averaged profiles of temperature and density 
from the three-dimensional simulated cluster data. We 
have therefore eliminated any systematic effect resulting 
from conversion of observed quantities. So we expect 
that the resulting "hydrostatic masses" should be closer 
to the true values than ones which would be observation- 
ally derived. We then compared these hydrostatic masses 
to the true mass for each cluster. We did this for a series 
of redshift intervals between 0.0 and 1.5, and separated 
clusters between CC and NCC. The result is shown in 
Figure [THI 

The average estimated cluster masses assuming hydro- 
static equilibrium are biased low for all the clusters by 
«15%. This bias is constant for CC clusters at differ- 
ent redshifts but appears to be slightly worse for NCC 
clusters at earlier epochs (~19%). In addition to the 
bias, the scatter in these mass estimates is high. 20% to 
30% underestimates are possible at the la level. Impor- 
tantly, and what is new here, CC clusters are no better 
than NCC clusters as biased mass indicators. Both are 
equally low, although the scatter in CC clusters is about 
h alf that of the NCC clusters. 

iMarkevitch fc Vikhlininl (|2007l ) similarly cast doubt on 
hydrostatic equilibrium in CC clusters due to the com- 
mon presence of observed cold fronts and the inferred 
"gas sloshing". For the clusters in our simulations, we 
find that the kinetic energy of bulk gas motions con- 
tributes at the »10% level compared to the total energy 
(see also lRasia et al.ll2006|) . 

If our numerical clusters are representative of real clus- 
ters, the apparent significant deviations from hydrostatic 
equilibrium for both CC and NCC clusters must be con- 
sidered in choosing to use them for precision cosmology 
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Fig. 13. — A sample of X-ray hardness ratio (HR) maps from 4 different simulated clusters (2 CC on top and 2 NCC on bottom) with 
fields of view of 0.33r2oo- The two circles in the upper left image delineate the "transition region" between 0.05r2oo to 0.3r2oo as described 
in Figure[7] Each HR map was made by dividing an image of the X-ray surface brightness from 2 to 8 keV by an image from 0.5 to 2.0 keV 
(i.e., typical Chandra hard and soft bands). Each is then normalized by the hardness ratio corresponding to the cluster virial temperature 
so that values < 1 represent gas with T < T viria i. Note that non-central cool "blobs" are generally infalling halos. 

estimators. 

6. SUMMARY AND CONCLUSIONS 

Galaxy clusters are complicated, generally non- 
equilibrium systems where nongravitational physics is 
important in the cores. To the best of our knowledge, no 
previous numerical simulations have been able to produce 
both cool core (CC) and non-cool core (CC) clusters in 
the same numerical volume. Our heating and cooling pre- 
scription (with an approximate balance between heating 
and cooling), however, has resulted in a simulation with 
both CC and NCC clusters. The temperature profiles of 
our numerical CC clusters qualitatively match observa- 
tions, although the central gas densities are higher than 
observed. Our fraction of cool cores is low compared to 
that of recent observed samples. On the other hand, the 
distinction in /3-model parameters (r c and 0) between 
CC and NCC clusters observed in samples of real clus- 
ters is reflected in our numerical clusters. Similarly, the 
distributions of gas fraction and gas mass with emission- 
weighted projected temperatures agree fairly well with 
observations. Overall, our numerical clusters have gen- 
eral characteristics that concur with X-ray data of ob- 




0.4 0.6 0.8 1.0 1.2 
Hardness Ratio (2-8 keV/0.5-2 keV) 

Fig. 14. — The cumulative fraction of pixels below each har dness 
ratio value for the sample of simulated clusters in Figure M3l The 
plot was assembled for regions between 0.05 — 0.3r"2oo (excludes cool 
core in CC clusters). The NCC clusters are centered roughly at 1, 
meaning that most of their gas is approximately equal to T v i r i a i. 
The CC clusters are noticeably cooler in this transition region. 
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Fig. 15. — Histograms of the ratio of the real space densities of 
halos surrounding numerical clusters with M200 = 1 — 6 X 1O 14 M0 
for CC versus NCC clusters. 
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Fig. 16. — The fractional deviation of cluster masses estimated 
assuming hydrostatic equilibrium versus the true virial masses as a 
function of redshift. The standard deviations in the distributions 
are also shown. 

served clusters. 

We propose an answer to the question posed in the ti- 
tle of this paper, i.e., only some clusters have cool cores 
because of evolutionary differences driven by early major 
mergers. Our numerical simulations suggest that the his- 
tories of cool core and non-cool core clusters are signifi- 
cantly different. Our NCC numerical clusters suffer early 
major mergers when nascent cool cores are destroyed. 
CC clusters, on the other hand, grow more slowly without 
early major mergers. CC clusters have a broad "transi- 
tion region" in their gas distribution extending between 



the cool core and a radius of ss 0.3r2oo where the gas 
fraction is higher than for NCC clusters and the temper- 
ature profile is nearly isothermal. This transition region 
and difference in evolution lead to a number of testable 
predictions for X-ray observations of real clusters. 

We find that the fraction of cool core clusters is a 
strong function of mass with fewer CC clusters at higher 
gas masses. This general trend qualitatively agrees with 
the analysis of recen t nea rby X-ray c l uster samples by 
lO'Hara et all (pOOlft and I Chen et all {2007ft . On the 
other hand, we do not find any significant variation in 
the fraction of numerical co ol cores with redshift in con- 
trast to the recent claim by Vikhlini n et al.l (|2006al) . 

The X-ray surface brightness profiles for NCC clusters 
are well fit by single /3-models whereas the outer emis- 
sion for CC clusters is biased low compared to /3-models. 
The resulting gas densities and gas masses of CC clusters 
estimated from single /?-model extrapolations are biased 
high by factors of 3-4. 

CC clusters have ~40% more cool gas beyond the cores 
within the transition region than do NCC clusters. This 
results in a very different distribution of X-ray hardness 
ratios beyond the cool core for CC versus NCC clusters. 
We predict that such differences will be observable with 
current X-ray imagers. 

There are some indications that the supercluster en- 
virons for CC and NCC clusters are different from each 
other today and in past epochs. At z > 1, NCC clusters 
appear to have more halos in their neighborhoods than 
CC clusters. At z < 0.3, this trend is reversed with more 
halos around CC clusters. This separation between CC 
and NCC clusters for low z clusters qualitatively agrees 
with supercluster density calculations for nearby Abell 
clusters. 

Finally, we find that both CC and NCC clusters are 
biased low in their mass estimation by ssl5% assuming 
hydrostatic equilibrium. In this sense, it appears that 
CC clusters are no better than NCC clusters as mass 
estimators, unlike what is generally assumed. This is 
important to consider in using CC clusters for precision 
estimations of cosmologic al parameters. 

In an upcoming paper ([Gantner et al.l [20071 ). we will 
compare the above predictions with X-ray observations 
of rich clusters from both the Chandra and ROSAT 
archives. The initial agreement is quite good. 

There are some important remaining issues with the 
current simulations. We plan to address these with a se- 
ries of new numerical simulations at higher resolution to 
overcome the limitations of the current computational set 
of clusters. We will explore the impact of Qb and eg on 
the fraction of cool cores that are produced in the compu- 
tational volume. We will also refine our heating/cooling 
prescription to better match current observational con- 
straints. Once the influence of these factors is understood 
on the creation of CC and NCC clusters, the fraction of 
cool core clusters could serve as an important new con- 
straint on cluster baryonic physics and/or dark energy 
models. 
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APPENDIX 

APPENDIX: X-RAY PROFILE FOR AN ADIABATIC ICM IN AN NFW CLUSTER POTENTIAL 

What does a good fit to the X-ray surface brightness profile by a /3-model impl y about the dynamic al state of the 
cluster gas when it is nonisothermal? When the intracluster gas is relaxed in a I Navarro etafl (Il99l (NFW) dark 
matter potential (derived from N-body simulations ) and in hydrostatic eq uilibrium, we can solve a simple equation for 
its radial distribution (see also lNavarro et aT]ll997l : [Makino fc Asanolll999l ). The equation for hydrostatic equilibrium 
is a simplification of the Euler equations for an ideal fluid, setting fluid velocity to zero gives 



gas 



-Pgas9, 



(Al) 



where P gas is the pressure of the ICM gas, p gas indicates the gas density, and g is the local gravitational acceleration. 
Under the assumption of spherical symmetry, we can simplify this to 



dP n , 



dr 



-Pgas{r)g{r). 



(A2) 



We assume here that the dark matter potential dominates the gravitation, and do not include the contribution of the 
gas to the potential, which should result in only minor error. In that case, we can write g(r) from an NFW dark 
matter profile as 

GM <R 

9{r) = (A3) 

where G represents the universal gravitational constant, and M <R indicates the dark matter mass inside the radius of 
interest. That mass can be calculated by integrating the NFW profile 



1 



where 



r 



(A4) 



(A5) 



and r c is the core radius and po.dm is the central normalization of the profile. Integrating the profile to get the total 
enclosed mass 

1 



M <R = 47rr c p , dri 



rX 2 dx 



x(l + x) 2 

We assume that the gas follows a nonisothermal, adiabatic equation of state such that 

P = kp 1 . 

where k is a constant. Then, the hydrostatic equilibrium equation to be solved can be written as 



(iP 

Uj± gas 



dr 



= -P 



l/ 7 
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POMn 
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x{\ 



l dx 



(A6) 



(A7) 



(A8) 



Numerical integration of Eq. 8 results in the profiles shown in Figure [TT] The solid line is the solution to Eq. 8, 
the dashed is the NFW dark matter density profile, and the dotted is a standard /3-model fit to the gas density. Note 
that the /3-model fits Eq. 8 very well for r > 0.2r c but the /3-profile is somewhat flatter in slope within the core in 
comparison to Eq. 8. This suggests that good /3-model fits to X-ray profiles imply gas that is relatively relaxed within 
the dark matter potential. 
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